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Abstract 

Perturbative QCD uses the Faddeev-Popov gauge-fixing procedure, which leads 
to ghosts and the local BRST invariance of the gauge-fixed perturbative QCD action. 
In the asymptotic regime, where perturbative QCD is relevant, Gribov copies can be 
neglected. In the nonperturbative regime, one must adopt either a nonlocal Gribov- 
copy free gauge (e.g., Laplacian gauge) or attempt to maintain local BRST invariance 
at the expense of admitting Gribov copies. These issues are explored. In addition, 
we discuss the relationship between recent Dyson-Schwinger based model calculations 
of the infrared behavior of QCD Green's functions and the lattice calculation of these 
quantities. 
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1. Introduction 



Perturbative quantum chromodynamics (QCD) is formulated using the Faddeev-Popov 
gauge-fixing procedure, which introduces ghost fields and leads to the local BRST invari- 
ance of the gauge-fixed perturbative QCD action. These perturbative gauge fixing schemes 
include, e.g., the standard choices of covariant. Coulomb and axial gauge fixing. These are 
entirely adequate for the purpose of studying perturbative QCD, however, they fail in the 
nonperturbative regime due to the presence of Gribov copies. Perturbative QCD works be- 
cause in doing a weak-field expansion around the = configuration these Gribov copies 
are not encountered.'^ 

One could define nonperturbative QCD by imposing a non-local Gribov-copy free gauge 
fixing (such as Laplacian gauge) or, alternatively, one could attempt to maintain local BRST 
invariance at the cost of admitting Gribov copies. One of the well-known difficulties for the 
latter option is the problem of pairs of Gribov copies with opposite sign giving a vanishing 
path integral.EJ'ISJ^J^ Whether or not a local BRST invariance for QCD can be maintained 
in the nonperturbative regime remains an open problem. 

The standard lattice definition of QCD is equivalent to the choice of a Gribov copy free 
gauge-fixing. There is a negligible chance of selecting two gauge-equivalent configurations 
(strictly zero except for numerical round-off error). Calculations of physical observables are 
unaffected by arbitrary gauge transformations on the configurations in the ideal gauge-fixed 
ensemble. A lattice QCD calculation using an ideal gauge-fixed ensemble will give a result 
for a gauge-invariant (i.e., physical) quantity which is identical to doing no gauge fixing at 
all, i.e., equivalent to the standard lattice calculation of physical quantities. 

We begin by reviewing the standard arguments for constructing QCD perturbation the- 
ory, which use the Faddeev-Popov gauge fixing procedure to construct the perturbative 
QCD gauge-fixed Lagrangian density. The naive Lagrangian density of QCD is Cqcb = 
— \F'^^F^y _l_ J2f 'lf{'^P~ '^f)^fi where the index / corresponds to the quark flavours. The 
naive Lagrangian is neither gauge-fixed nor renormalized, however it is invariant under local 
SU{'i)c gauge transformations g{x). For arbitrary, small uj°'{x) we have 
g[x) = exp {—igs {X"" /2) uj''{x)} G SU{3), where the A'*/2 = t°- are the generators of the 
gauge group SU{3) and the index a runs over the eight generator labels a = 1, 2, 8. 

Consider some gauge-invariant Green's function (for the time being we shall concern 
ourselves only with gluons) (J7| T{d[A]) \Q) = JVA 0[A] e^^t^V / VA e^^t^l , where 0[A] 
is some gauge- independent quantity depending on the gauge field, A^{x). We see that the 
gauge- independence of 0[A] and S[A] gives rise to an infinite quantity in both the numerator 
and denominator, which must be eliminated by gauge-fixing. The Minkowski-space Green's 
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functions are defined as the Wick-rotated versions of the Euchdean ones. 

The gauge orbit for some configuration is defined to be the set of all of its gauge- 
equivalent configurations. Each point Af^ on the gauge orbit is obtained by acting upon 
with the gauge transformation g. By definition the action, S[A], is gauge invariant and so 
all configurations on the gauge orbit have the same action, e.g., see the illustration in Fig. ^ 



gauge A.^ 
orbit 

Fig. 1. Illustration of the gauge orbit containing and indicating the effect of acting on A^ with 
the gauge transformation g. The action S[A] is constant around the orbit. 



Any gauge-fixing procedure defines a surface in gauge-field configuration space. Fig. |21 is 
a depiction of these surfaces represented as dashed lines intersecting the gauge orbits within 
this configuration space. Of course, in general, the gauge orbits are hypersurfaces as are 
the gauge-fixing surfaces. Any gauge- fixing surface must, by definition, only intersect the 
gauge orbits at distinct isolated points in field configuration space. For this reason, it is 
sufficient to use lines for the simple illustration of the concepts here. An ideal (or complete) 
gauge-fixing condition, F[A] = 0, defines a surface called the Fundamental Modular Region 
(FMR) that intersects each gauge orbit once and only once and typically where possible 
contains the trivial configuration = 0. A non-ideal gauge-fixing condition, F'[A] = 0, 
defines a surface or surfaces which intersect the gauge orbit more than once. These multiple 
intersections of the non-ideal gauge fixing surface(s) with the gauge orbit are referred to as 
Gribov copies .'^''^''^''^ Lorentz gauge = 0) for example, has many Gribov copies 

per gauge orbit. By definition an ideal gauge fixing is free from Gribov copies. The ideal 
gauge- fixing surface F[A] = specifies the FMR for that gauge choice. Typically the gauge 
fixing condition depends on a space-time coordinate, (e.g., Lorentz gauge, axial gauge, etc.), 
and so we write the gauge fixing condition more generally as -F([v4]; x) = 0. 

Let us denote one arbitrary gauge configuration per gauge orbit, ^4° , as the origin for that 
gauge orbit, i.e., corresponding to (7 = on that orbit. Then each gauge orbit can be labelled 
by A° and the set of all such A^ is equivalent to one particular, complete specification of 
the gauge. Under a gauge transformation, g, we move from the origin of the gauge orbit 




§2. Gribov Copies and the Faddeev-Popov Determinant 



3 



A,=0 

Fig. 2. Ideal, F[A], and non-ideal, -F'[^], gauge-fixing. 



to the configuration, ^4^, where by definition A° — ^ — gA^^g^ — {i/ gs){dfj_g)gK Let us 
denote for each gauge orbit the gauge transformation, g = g[A'^], as the transformation 
which takes us from the origin of that orbit, to the corresponding configuration on the 
FMR, A™^ = A^, which is specified by the ideal gauge fixing condition F([A^];a;) = 0. In 
other words, an ideal gauge fixing has a unique g which satisfies x)\g = and hence 

specifies the FMR as A^ = A™^ e FMR. Note then that we have JVA^ J VA° JVg ^ 
jVA''^'^jV{g-~g). 

The inverse Faddeev-Popov determinant is defined as the integral over the gauge group 
of the gauge-fixing condition, i.e., 

-1 



A-/[A^^^] ^Jvg S[F[A]] = J Vg 5{g - ~g) |det (^^^^ 



(2-1) 



Let us define the matrix Mf[A] as MF{[A];x,yY'' = SF''{[A];x)/Sg\y) . Then the Faddeev- 
Popov determinant ioT an arbitrary configuration Af^ can be defined as = |det Mi?[74]|. 
(The reason for the name is now clear) . Note that we have consistency, since Aj} [A™^] = 
A-'[A^]^JVg S{g - g)A-/[A]. 

We have 1^ fVg Af[A] S[F[A]] ^ fV{g- g) Af[A\ 5[F[A\\ by definition and hence 

j PA™^= j VA™""! V{g - g) Af[A]S[F[A]] VA Af[A]5[F[A]] (2-2) 

Since for an ideal gauge-fixing there is one and only one g per gauge orbit, such that 
F([74];a;)|g = 0, then |detMi?[yl]| is non-zero on the FMR. It follows that since there is 
at least one smooth path between any two configurations in the FMR and since the determi- 
nant cannot be zero on the FMR, then it cannot change sign on the FMR. The first Gribov 
horizon is defined to be those configurations with detMj7'[yl] = which lie closest to the 
FMR. By definition the determinant can change sign on or outside this horizon. Clearly, 
the FMR is contained within the first Gribov horizon and for an ideal gauge fixing, since 
the sign of the determinant cannot change, we can replace | detMi;'] with detM^T^, [i.e., the 
overall sign of the functional integral is normalized away in the ratio of functional integrals] . 

These results are generalizations of results from ordinary calculus, where 
|det {df,/dx,)\-i^ = jdx,--- dxn S^-\f{x)) 
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and if there is one and only one x which is a solution of f{x) =0 then the matrix Mij = 
dfi/dxj is invertible (i.e., non-singular) on the hypersurface f{x) = and hence detM 7^ 0. 

§3. Generalized Faddeev- Popov Technique 

Let us now assume that we have a family of ideal gauge fixings F([74];a;) = — 
c{x) for any Lorentz scalar c{x) and for being some Lorentz scalar function, (e.g., 

d^Aij,{x) or n^Aij,{x) or similar or any nonlocal generahzations of these). Therefore, using 
the fact that we remain in the FMR and can drop the modulus on the determinant, we have 
J VA™^ = / VA det Mf[A] S[f[A] - c] . Since c{x) is an arbitrary function, we can define 
a new "gauge" as the Gaussian weighted average over c{x), i.e., 

J VA™"" <x Jvc exp y d'^xcixf^ J VA detM^f^l] S[f[A] - c] 

(X JvA detMF[A]expi^-^J d^xf{[A];x)^^ 

xe:<pl-^j d*xf([A];xA , (3-1) 

where we have introduced the anti-commuting ghost fields x X- Note that this kind 
of ideal gauge fixing does not choose just one gauge configuration on the gauge orbit, but 
rather is some Gaussian weighted average over gauge fields on the gauge orbit. We then 
obtain 

where 

SdQ, Q, A x,x]-J d'x -\f'^^'F;, - ^ (/([^]; x))^ + ^ qfiilp- mf)qf 

+ J d^xd^y x{x)MF{[A]-x,y)x{y) ■ (3-3) 

§4. Standard Gauge Fixing 

We can now recover standard gauge fixing schemes as special cases of this generalized 
form. First consider standard covariant gauge, which we obtain by taking /([74];a;) = 
dnA^^{x) and by neglecting the fact that this leads to Gribov copies. We need to evalu- 
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ate Mp[/1] in the vicinity of the gauge- fixing surface (specified by g): 
MF{[A];x,y) 



^ SF%[A];x) ^ 5[d,A->^{x)-c{x)] ^ 

8g\y) 5g\y) ^ 5g\y)] ' 



(4-1) 



Under an infinitesimal gauge transformation about the FMR, 6g = g — g, we have {A^)ij, 
where 



{A^^'%{x) = {A%{x) + g,r"^u;\x)Al{x) - d,u;^{x) + 0{ 



(4-2) 



and hence near the gauge fixing surface (i.e., for smaU fluctations along the orbit around 
A™^) using MF{\A]-x,yY^ = d^[6A^''{x)/6{6uj\y)])\^=o we find 

MFi[A];x, yr' = ( [-d^'^S'^' + g,r'''A'^^ix)]5^'\x - y)) . 

We then recover the standard covariant gauge-fixed form of the QCD action 

Sdq, q, A, X,X]-J d'x —F'^^^'F;, - 1 {d.A^f + ^ qf{ip- mf)qf 

Hd^XaWS'^'-gfabcA^^Xb. (4-3) 

However, this gauge fixing has not removed the Gribov copies and so the formal manipu- 
lations which lead to this action are not valid. This Lorentz covariant set of naive gauges 
corresponds to a Gaussian weighted average over generalized Lorentz gauges, where the 
gauge parameter ^ is the width of the Gaussian distribution over the configurations on the 
gauge orbit. Setting ^ = we see that the width vanishes and we obtain Landau gauge 
(equivalent to Lorentz gauge, d^An{x) — 0). Choosing ^ = 1 is referred to as "Feynman 
gauge" and so on. We can similarly derive the QCD action for axial gauge. 

We can similarly recover the standard QCD action for the axial gauges, where n^A'*(x) = 
0. Proceeding as for the generahzed covariant gauge, we first identify /([>!]; x) = nij,A^{x) 
and obtain the gauge-fixed action 



SdQ.q,A] 



(4-4) 



Taking the "Landau-hke" zero- width hmit ^ ^ we select n^A^^{x) = exactly and recover 
the usual axial-gauge fixed QCD action. Axial gauge does not involve ghost fields, since in 
this case 



M^([>l™^];a;,y)«''=n^ 



ab 



5A''>'ix) 



8uj\y)] 



-n^{[-d^''5^'']5^^\x-y)), 



(4-5) 



uj=0 
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which is independent of since n^y4™^(x) = 0. In other words, the gauge field does 
not appear in Mi7[yl] on the gauge- fixed surface. Unfortunately axial gauge suffers from 
singularities which lead to significant difficulties when trying to define perturbation theory 
beyond one loop. A related feature is that axial gauge is not a complete gauge fixing 
prescription. While there are complete versions of axial gauge on the periodic lattice, these 
always involve a nonlocal element, or reintroduce Gribov copies at the boundary so as not 
to destroy the Polyakov loops in the axial gauge-fixing direction. 

§5. Discussion and Conclusions 

There is no known Gribov-copy-free gauge fixing which is a local function of In 
other words, such a gauge fixing cannot be expressed as a function of and a finite 

number of its derivatives, i.e., F([y4];a;) 7^ F{d^, A^{x)) for all x. Hence, the ideal gauge- 
fixed action, S*^ [■■■], in Eq. (j3-3j) becomes non-local and gives rise to a nonlocal quantum field 
theory. Since this action serves as the basis for the proof of the renormalizability of QCD, the 
proof of asymptotic freedom, local BRST invariance, and the Dyson-Schwinger equation^''^ 
(to name but a few) the nonlocality of the action leaves us without a first-principles proof 
of these features of QCD in the nonperturbative context. 

The lattice implementation of Landau gauge finds local minima of the gauge fixing func- 
tional, which correspond to configurations lying inside the first Gribov horizon. The re- 
maining Gribov copies after this partial gauge fixing then necessarily all have the same sign 
(positive) for the Faddeev-Popov determinant and hence add coherently in the functional 
integral. This ensures that the ghost propagator is positive definite.^''^ The derivation of 
the Dyson-Schwinger equations is based on the fact that the integral of a total derivative 
vanished provided that the surface integral of the integrand vanishes when integrated over 
the boundary of the region. Since the Faddeev-Popov determinant vanishes on the first 
Gribov horizon, then we can still derive the standard Dyson-Schwinger equations from the 
Landau gauge fixed QCD action even if we restrict the gauge fields to lie within the first 
Gribov horizon. This is equivalent to requiring that the ghost propagator be positive defi- 
nite. Thus it is valid to compare lattice Landau-gauge calculations with Dyson-Schwinger 
based calculations (with a positive definite ghost propagator), since these both consist of 
considering configurations within the first Gribov horizon. An extensive body of lattice cal- 
culations exist for the Landau gauge gluorP and quar)^^ propagators and most recently 
for the quark-gluon vertex^ Similarly, calculations in Laplacian gauge (an ideal gauge) 
fixing have also recently become available .'^''^ 

It is well-established that QCD is asymptotically free, i.e., it is weak-coupling at large 
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momenta. In the weak coupling limit the functional integral is dominated by small action 
configurations. As a consequence, momentum-space Green's functions at large momenta will 
receive their dominant contributions in the path integral from configurations near the trivial 
gauge orbit, i.e., the orbit containing = 0, since this orbit minimizes the action. If we use 
standard lattice gauge fixing, which neglects the fact that Gribov copies are present, then at 
large momenta J T>A will be dominated by configurations lying on the gauge-fixed surfaces 
in the neighbourhood of each of the Gribov copies on the trivial orbit. Since for small 
field fiuctuations the Gribov copies cannot be aware of each other, we merely overcount the 
contribution by a factor equal to the number of copies on the trivial orbit. This overcounting 
is normalized away in the ratio of functional integrals. Thus it is possible to understand why 
Gribov copies can be neglected at large momenta and why it is sufficient to use standard 
gauge fixing schemes as the basis for calculations in perturbative QCD. 
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